Analysis Of Measured Transport Properties Of 
Domain Walls In Magnetic Nanowires And Films 



L. BERGER 

Physics Department, Carnegie Mellon University, Pittsburgh, PA 

15213 



Abstract: Existing data for soft magnetic materials of critical current for 
domain-wall motion, wall speed driven by a magnetic field, and wall electrical 
resistance, show that all three observable properties are related through a single 
parameter: the wall mobility fi. The reciprocal of fi represents the strength 
of viscous friction between domain wall and conduction-electron gas. And /i is 
a function of the wall width, which depends in turn on the aspect ratio t/w, 
where t and w are the thickness and width of the sample. Over four orders of 
magnitude of /z, the data for nanowires show fx oc (t/w)~ 2 - 2 . This dependence is 
in approximate agreement with the prediction of the 1984 Berger theory based 
on s-d exchange. On the other hand, it is inconsistent with the prediction of 
the 2004 Tatara and Kohno theory, and of the 2004 Zhang and Li theory. 

PACS numbers: 73.50.Jt, 75.60.Ch, 75.70.-i; Keywords: domain wall, mag- 
netic nanowire, galvanomagnetic phenomena. 

I. INTRODUCTION 

Many years ago, we predicted 1 that a dc current would apply a force on 
a domain wall, through the s-d exchange interaction. We also observed 2 ' 3 
wall displacements caused by such a force, in Ni-Fe films containing Neel walls. 
As predicted, the wall displacements had the same sense as the drift speed 
of the electron-like carriers. Can et al. 4 saw similar current-induced wall 
displacements in Ni-Fe films with smaller patterned width of 20 \xm. Since 
then, the same kind of results have been reported by many authors for magnetic 
nanowires 5 ~ 14 . These nanowires, mostly made of Ni-Fe, have a much larger 
ratio of thickness t to width w, resulting in thinner walls and, consequently , 
to a stronger interaction of wall and electron gas. The values of t and w, of 
the measured critical current density j c where wall motion starts in zero field, 
and of the measured coercive field H c for wall motion at zero current are the 
important measured quantities. In some cases, j c or H c have to be inferred from 
more indirect information given in the paper. 

Ebels et al. 15 as well as Khizroev et al. 16 and Lepedatu and Xu 5 have 
measured the extra electrical resistance Rdw caused by a wall, in nanowires. 
The Ebels cobalt sample had a circular cross section with t=w as the diameter. 

The wall speed v w induced by an easy-axis field H is proportional to H, and 
the ratio \x = v w /H is called the wall mobility. Patton and Humphrey 17 as well 
as Shishkov 18 have measured [i for walls in Ni-Fe films over a range of thickness 
t. They also summarized earlier work. We use their results for t < lfim, where 



eddy-current damping is negligible. Freitas and Berger 2 have measured /i in 
Ni-Fe in the same samples where they moved walls with a current. Finally, Ono 
et al. 19 measured fj, on a Fe-Ni nanowire. 

The purpose of the present paper is to show that the concept of wall mobility 
can be used to characterize in a consistent manner the results of experiments 
of current-induced wall motion, or of measurements of wall resistance, just as 
well as those of field-induced wall motion. No matter the origin of a value of 
H, we find that value to be the same function of the ratio t/w of thickness 
to width of the sample. The physical reason 1 is that in all three cases we 
have viscous friction between a domain wall and the electron gas. The value 
of the friction force between the two is the same whether the electron gas is 
moving with respect to an initially static wall (experiments of wall resistance 
or of current-induced wall motion) , or whether the wall moves and the gas is at 
rest (field- induced wall motion). 

We also show that these experiments agree with the 1984 Berger theory, but 
not with the 2004 theories by Tatara and Kohno or by Zhang and Li. 

II. WALL MOBILITY DERIVED FROM 
CRITICAL CURRENT OR WALL RESISTANCE 



We start with the force Fh exerted by an easy-axis field H on a 180° wall, 
given 20 by 

F H = 2fj, M s Htw. (1) 

In the steady state, this force must be balanced by a viscous damping force 
Fd, giving Fh + Fd = 0. Combining this equation and Eq. (1) with the mobility 
definition \x = v w /H given earlier, we obtain 

F d = -2pt M s tw^v w . (2) 

In the 1984 Berger theory 1 , the origin of this damping force is the friction 
of the moving wall with the electron gas at rest. If, instead, the gas is moving 
because an electric current density j is present, the average drift speed of the gas 
is v e = —j/n e e = jR , where n e is the electron density in a one-band model, 
and Rq the ordinary Hall coefficient. If the interaction depends only on relative 
motion, we expect 1 v e to generate a drive force F e on the wall similar to the 
damping force Fd generated by v w (see Eq. (2)), though of the opposite sign: 

F e = 2 f i Q M s tw f i' 1 (3 1 v e . (3) 

Here, (3\ is 1 a correction factor found experimentally 2,3 to be equal to about 
2 for Ni-Fe films. 

There is also a pinning force F p applied to the wall by lattice defects. By 
analogy with Eq. (1), it is 

F p = ±2/j, M s H c tw. (4) 

where H c is the coercive field. At the critical current density j c for wall 
motion, F e just balances F p , with F e + F p — 0. Combining this equation and 
the relation v e = Roj with Eqs. (3) and (4), we obtain finally: 
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From the measured j c ,H c obtained from Refs. 2-14, we find \x for films 
and nanowires by Eq. (5), assuming 2 Rq = —1.36 x 10~ 10 ra 3 /C for Ni 80 Fe2o 
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and 21 -1.1 x l(r 10 m 3 /C for Co, also 21 -1.25 x lir^m^/C for Fe 5Q Co 50 . 
We plot /i versus t/w in Fig. 1 for all these samples. In several cases, the 
quantities actually measured were not j c , H c but critical-current variations Aj c 
corresponding to field variations AH. But Eq. (5) keeps the same form. 

To the drive force F e (Eq. (3)) exerted by electrons on the wall corresponds 
a reaction force — F e applied to the electron gas, which manifests itself 22 as a 
wall resistance Rdw > 0: 

2(3 1 ^M s Rl 

Rdw = • (O) 

fiwt 

From the measured Rdw for nanowires 5 ' 15 ' 16 , we find \x by Eq. (6), and plot 
it versus t/w in Fig. 1. We assume /j.qM s = IT for NisoFe 2 o and 1.8 T for Co, 
and the same 0i,Ro as before. Note that Eq. (6) does not hold in the case of 
Rdw mechanisms where the momentum flows from the electrons to the lattice, 
rather than to the wall itself. One example is the mechanism of Gregg and of 
Levy and Zhang 23 . Another is the anisotropic magnetoresistance. Because of 
this, we have only used Rdw measurements where the Rdw value predicted by 
Eq. (6) is larger than the Gregg or anisotropic-magnetoresistance ones, i.e., for 
small A corresponding to large t/w > 0.1. Indeed, the Gregg resistance 23 
varies like Ag 1 oc {t/w) 1 / 2 , slower than our Rd w of Eq. (6). 

III. EMPIRICAL RELATION BETWEEN 
WALL MOBILITY AND ASPECT RATIO 

The wall mobilities \i obtained by the three different methods are plotted 
in Fig. 1 versus t/w, for all films and nanowires, on double logarithmic scales. 
The fact that the directly measured wall mobilities (one case for nanowires and 
many for films) are roughly consistent with those derived from critical current 
(many for nanowires and three for films) and with those from wall resistance 
(six for nanowires) indicates that the phenomenological idea 1 of viscous friction 
between electron gas and domain wall, used in the preceding data analysis, is 
sound. 

We see (Fig. 1) that the \x values for nanowires start at the same level as 
those for films, but plunge with a large negative slope, dropping by four orders 
of magnitude as t/w increases. Although the dispersion of the data points is 
rather large, the best fit of a straight line to the nanowire data is with a slope 
of -2.2 (dashed line in Fig. 1). Only the data of Klaui et al. 14 really differ from 
this scheme. The average horizontal deviation of data points from the dashed 
line is 2.8 decibels, excluding the point of Ref. 14. This compares to a total 
range of variation of t/w by a factor of 100, or 20 decibels, leading to an accuracy 
of about 14% on the slope of the dashed line. The correlation coefficient of the 
nanowire data is found to be 80%. 

The film data (Fig. 1), including those of Gan et al. (Ref. 4), seem to lie 
in a nearly horizontal band. The reason is probably that t/w is too small for 
the shape anisotropy to dominate over the induced anisotropy, so that the wall 
width varies little with t/w, except near t/w = 10~ 5 where the walls change 
from Nccl to Bloch type. 

We still need to find the relation between Ao and t /w in the case of nanowires, 
where the main contribution to the anisotropy energy is 25 the shape anisotropy. 
The latter is connected with magnetic poles at the wire surface. For the sim- 
plest kind of head-to-head 180° domain wall with no vortices, we obtain A by 
minimizing the sum of exchange and anisotropy energies: 
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A ° (|arcton(^)) 1 /2- (7) 

Here, l ex — (2A/ hqM^) ~ 6nm is the exchange length of Ni$oFe2o, and 
A ~ 1.5 x lO~ 11 J/m the exchange stiffness. The wall mobility is often written 
24 in terms of the Gilbert damping parameter a as fi — jAo/a, where 7 is the 
gyromagnctic factor. The results of Fig. 1 for nanowires imply that, like /i, a 
is dependent on t/w through the wall width Ao, with a oc (t/w) 1 - 7 oc Aq 3 ' 4 if 
t/w < 1. 

IV. MICROSCOPIC THEORY OF VISCOUS FORCES 

We now turn to the 1984 microscopic theory 1 , which was based on the 
force applied to conduction-electron spins by the gradient of the s-d exchange 
field present in the wall. This force between electrons and wall is similar to the 
one involved in the famous Stcrn-Gcrlach experiment, but with the magnetic 
field replaced by an exchange field. It is important not to confuse this force 
(momentum transfer) with the spin-transfer torque, which is not equivalent to 
a force and was discussed in the 1986 paper of Ref. 22. The reciprocal of the 
wall mobility is predicted to be: 

n r+00 c\ 
^ = 5^ / dx[Uv,{x) - V^x)?. (8) 



16/u M s J_ 00 dx 

Here, r-j-j, is 1 the relaxation time of electrons between bands. It is largest in 
the particular case of spin-up and spin-down bands, considered here, where it 
becomes the spin-relaxation time. One of the two d(V^(x) — Vi(x))/dx factors, 
when combined with the 774 factor, represents a non-equilibrium transfer of 
electron density between the two bands, i.e., a spin accumulation. This spin 
accumulation is a scalar which has 1 opposite signs on the two sides of the 
domain wall. Also, D^ s is the density of states of 4s conduction electrons, 
assumed the same for spin-up and spin-down states. And V-\(x) > 0, V[(x) < 
are potential-energy functions 1 of maximum order 10 _3 eV which vanish outside 
the wall. If spin-up and spin-down conduction electrons have the same Fermi 
velocity vf, we have 



V T (x) = -Vl(x) = J sd S(l-cosj(x));-f(x) = arctan ( 2SJ ^31/2 ^)- ( 9 ) 

Here, J s d is the s-d exchange integral, and S the atomic spin of magnetic 3d 
electrons assumed localized. Also, 6(x) is the angle of the local magnetization 
in the wall with the domain magnetization, given 20 by x/Ao = In tan(9/2). 
Finally, 7(x) is the local angle between conduction-electron spin and magnetic 
spin in the wall, later called "mistracking angle" by Gregg 23 . Semiclassical 
equations of motion 1 for conduction-electron spins are used to derive Eqs. (9). 
While the approximation 7(2:) <C lrad. was made in Ref 1, Eqs. (9) are valid 
even when 7(2;) ~ lrad.. Eqs. (9) show that V\ and arise from the existence 
of the mistracking angle. 

It was shown in Ref. 1 that the Stern-Gerlach forces on individual electrons 
of spin up and spin down can be written in the form of — dV^/dx, —dVi/dx, 
respectively. This leads to the second d(V](x) — Vi(x))/dx factor in Eq. (8). 
The spin accumulation mentioned above is needed if the total force created by 
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the left half of the wall is not to cancel the total force of the right half (see Fig. 
1 of Rcf. 1). 

We evaluate p numerically as a function of t/w by combining Eqs. (7-9). 
We use values J s d — 0.6eU, and poM s = IT, S = 0.5 for Ni &0 Fe2o- We derive 
the spin-relaxation time tjj. from the spin-diffusion length l sr = [2t-\J e 2 D{p^ + 
pi)) 1 ^ 2 where D is the full density of states including 3d electrons, and p\,P{ 
the spin-up, spin-down resistivities. The electronic specific heat 26 of Ni so Fe2o 
gives D = 1.09 xlO 48 J _1 m- 3 . We assume pi+pi = 133 x 10" 8 Om for Ni so Fe 2 o 
at room temperature, somewhat higher than the 105 x 10~ 8 f2m derived from the 
measured low-temperature p^,pi values 27 . Finally, using a value l sr — 60nm 
for a ferromagnet such as cobalt 28 , we obtain r-f j. = 6.7 x 10~ n s. 

If we use free-electron values vf = 1-57 x 10 e m/s, D^ s = 11.4 x 10 46 J -1 ™ -3 , 
suitable 29 for conduction electrons in Cu, then Eqs. (7-9) yield much too large 
p values for given t/w, corresponding to very small viscous coupling between 
electrons and wall. In order to fit the nanowire data of Fig. 1, we are forced to 
use Vf = 22.2 x 10 6 m/s, about 14 times larger than the free-electron value above. 
To be consistent, we have to reduce Di s by the same factor, toS.lxlO 45 ,/- 1 
The p predicted by Eqs. (7-9) are then plotted versus t/w as the solid curve in 
Fig. 1. The too large needed vf value mentioned above may represent a genuine 
band-structure effect in nanowires. But if, instead of l sr = 60nm, we used the 
surprisingly short ~ 4 — 5nm measured 27 in Ni-Fe films, an even larger vf value 
would be needed. 

We see (Fig. 1) that the fit between theory (solid curve) and the average 
of experiments (dashed line) is rather good. In the limit 7(2;) <C 1, realized at 
small t/w, Eq. (9) gives 7(2;) oc A^ 1 and V^Vj. oc Aq 2 . Plugging this into 
Eq. (8), we obtain p oc Aq. Using Eq. (7), this means p oc (t/w)~ 2 - 5 , giving a 
predicted slope is -2.5. As t/w and 7(2;) increase, the slope of the theoretical 
curve becomes less negative and is close, overall, to the average experimental 
slope of -2.2 indicated by the dashed line (Fig. 1). The mistracking angle 7(2:) 
changes from 0.19 rad. to 1.0 rad. as t/w increases over the range covered 
by the nanowires. Note that p oc (t/w)~ 2 - 2 implies p oc A 4 - 4 by Eq. (7), if 
t/w < 1. 

V. IMPORTANCE OF RANDOM SCATTERING 

In NisoFe2Q films or nanowires, the average of spin-up and spin-down elec- 
tron mean free paths has been estimated 27 at ~ 2nm. In Co, it is not much 
longer. This justifies the 1984 Bcrger theory , where the mean free path is 
assumed smaller than the wall width A = 8 — 70nm, given by Eq. (7) for the 
present nanowires. The small mean free path breaks up the coherence of the 
wavefunction, leading 1 to the validity of a semiclassical transport theory with 
a local conductivity defined at each point inside the wall 

Tatara and Kohno's recent theory 30 , like the Berger theory , is based on 
the "Stern-Gerlach force" caused by the gradient of the exchange field in the 
wall. However, they do not have any random electron scattering by solutes or 
phonons, so that their wavefunction is highly coherent over the whole wall. This 
leads to nearly-complete destructive interference between wavelets reflected at 
different points of the wall, resulting in a very small total force which varies 
exponentially 30 with wall width Ao, not as a power law as for the data of Fig. 
1 or for the 1984 Berger theory. This is shown as the dotted curve labelled T-K 
in Fig. 1, assuming J s d — 0.1eV,vp = 1-57 x 10 6 m/s. 

The theory by Thiaville et al. 31 is based on the same phenomenological 
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idea as our data analysis (Eqs. (2-3)), namely that the electron drift speed v e 
generates friction forces of the same magnitude as the wall speed v w . 

Zhang and Li 32 predict a force proportional to Aq" 1 . This is equivalent to 
a wall mobility [i oc A oc (t/w)~ - 5 , shown as a dotted line of slope -0.5 in 
Fig. 1, labelled Z-L. This is clearly inconsistent with the /i oc (t/w)~ 2 - 2 found 
experimentally by us in nanowires (dashed line in Fig. 1), and predicted by the 
Bcrger theory (solid curve in Fig. 1). In addition, the [i values predicted by 
Zhang and Li for values r-fj. = 0.7 x 10 _13 s and J s d = 0.6eV (Fig. 1) are too 
large by at least three orders of magnitude to explain the experimental results 
for t/w > 0.1. They correspond to a very small force on the wall. As in the 
case of the Tatara and Kohno theory, we chose the above parameter values in 
such a way as to maximize the force. 

Very recently, we learnt about a calculation by Zhao et al. 33 of the resistance 
of a domain wall at a point contact in the diffusive limit. Using our Eq. (6), 
their result is equivalent to a mobility \i oc A 2 , oc (t/w)^ 1 . 

VI. CONCLUSIONS 

We found that measured domain- wall electrical resistance and current-indu- 
ced as well as field-induced wall motion are all related through the concept of 
wall mobility jj,, and lead to the same dependence /i oc (i/u>)~ 2,2 on the aspect 
ratio t/w of the cross section of a magnetic nanowire (Fig. 1). Note that, for 
domain- wall resistance, this applies only at t/w > 0.1, where the resistance 
mechanism of Eq. (6) dominates. This dependence of /i on t/w agrees with 
the prediction (Eqs. (7-9)) of the 1984 Berger theory 1 , though a too large 
value of the Fermi velocity is needed to explain the observed /x values. But it is 
inconsistent with the 2004 Tatara and Kohno theory, as well as with the 2004 
Zhang and Li theory. 

I am grateful to Sonali Mukhcrjcc for useful discussions on this topic. 
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FIGURE CAPTIONS 

FIG. 1: Wall mobility li versus aspect ratio t/w of cross section of sample. 
Plain circles represent \i values from measurements of critical current density j c . 
Circles with a tail to the left represent li values from measurements of domain- 
wall resistance Rdw, and those with a tail to the right direct measurements of 
li. The numbers in the circles indicate the references from which the data come. 
The dashed straight line is fitted to the nanowire data, and has a slope of -2.2. 
The solid curve is the prediction of the 1984 Berger theory (see Eqs. (7-9)). 
The two dotted curves are the predictions of the Tatara and Kohno theory and 
of the Zhang and Li theory. 
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